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On pure states of n quantum bits, the concurrence entanglement monotone returns the norm of the inner product 
of a pure state with its spin-flip. The monotone vanishes for n odd, but for n even there is an explicit formula for 
its value on mixed states, i.e. a closed-form expression computes the minimum over all ensemble decompositions 
of a given density. For n even a matrix decomposition v = k\ aki of the unitary group is explicitly computable 
and allows for study of the monotone's dynamics. The side factors k\ and ki of this Concurrence Canonical 
Decomposition (CCD) are concurrence symmetries, so the dynamics reduce to consideration of the a factor. This 
unitary a is diagonal on a basis of GHZ-like states, with dynamics of the entire state space depending on the 
relative phases within a. In this work, we provide an explicit numerical algorithm computing v = k\ aki for n 
odd. Further, in the odd case we lift the monotone to a two-argument function. The concurrence capacity of v 
according to the double argument lift may be nontrivial for n odd and reduces to the usual concurrence capacity in 
the literature for n even. The generalization may also be studied using the CCD, leading again to maximal capacity 
for most unitaries. The capacity of v S3 h is at least that of v, so odd-qubit capacities have implications for even- 
qubit entanglement. The generalizations require considering the spin-flip as a time reversal symmetry operator in 
Wigner's axiomatization, and the original Lie algebra homomorphism defining the CCD may be restated entirely 
in terms of this time reversal. The polar decomposition related to the CCD then writes any unitary evolution 
as the product of a time-symmetric and time-antisymmetric evolution w.r.t. the spin-flip. En route we observe 
a Kramers' nondegeneracy: the existence of a nondegenerate eigenstate of any time reversal symmetric ji-qubit 
Hamiltonian demands (i) n even and (ii) maximal concurrence of said eigenstate. We provide examples of how 
to apply this work to study the kinematics and dynamics of entanglement in spin chain Hamiltonians. 
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I. INTRODUCTION 

The entanglement theory of two quantum bits is now well understood. For let p be a mixed two-qubit quantum 
state, described by a 4 x 4 Hermitian density matrix. Hill and Wootters (H97) describe all classes of p up to 
evolution by unitaries in terms of the concurrence. This concurrence is explicitly a function of the eigenvalues of 
p(o y )® 2 p(G y )® 2 , where the factor p = (o > )® 2 p(a > )® 2 may be interpreted as the spin-flip of p. Further, for pure 
states the entropy of the partial trace to either one-qubit subsystem is a one-to-one function of the concurrence, 
so that both measures agree as to which two-qubit states are more or less entangled. Local states (tensors) are 
unentangled, while states locally equivalent to Bell states have maximal entropy and concurrence. 

For other systems, entanglement theory is more complicated. Even for two d-level systems (qudits,) it is 
not typical to use a single function to quantify entanglement (V00), and research into generalized concurrences 
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continues (G04). Instead we focus on the multi-partite qubit case. The key point is that it is not sensible in n-qubits 
to speak of a unique maximally entangled state. More precisely, suppose now p is a 2" x 2" Hermitian density 
matrix describing a mixed n-qubit state. A unitary evolution is given by a 2" x 2" unitary matrix, say v, with the 
evolution being p i— > vpv^. A partial ordering of such p as more or less entangled follows by stipulating that (i) for 
v = ®" = i v j local unitary, p and vpv^ are equally entangled, while (ii) the p becomes no more entangled on average 
after applying any sequence of local measurements and local unitaries, i.e. after applying local completely-positive 
maps. More entangled is a partial order which has distinct maximal elements for n > 3. For example, in three 
qubits, two states which are maximally entangled yet locally inequivalent are given as follows (DVC01): 

|GHZ) = (1/V2)[|000) + |111)], |W) = (1/V3)[|001) + |010} + |100}] (1) 

There are nine distinct maxima of the partial order in four qubits (VDMD02), and strong theoretical evidence 
suggests that the number of such entanglement types grows quite rapidly with n (e.g. (MW02b).) 

To quantify multi-partite entanglement, one often uses functions known as entanglement monotones (V00; 
BL01). All such monotones must vanish on any local state. A monotone might also vanish on certain entangled 
states but definitively reports that a state is not local should its value be nonzero. The value on a mixed state p is 
defined to be the minimum over all ensemble decompositions of p of the ensemble weighted-average. A monotone 
is convex on density matrices, since entanglement does not increase under mixing of states. Monotones are also 
nonincreasing on average under local quantum operations and classical communication. Among popular mono- 
tones are Meyer's g-measure (MW02a; B03), the Schmidt measure (EB01), and certain polynomial invariants 
(BL01) of eigenvalues of density matrices representing stochastic mixtures of pure data states. 

The n-qubit concurrence is an entanglement monotone. To define the monotone, we first note that throughout 
13 refers to the spin-flip of the n-qubit state space. Concurrence for a pure state (WC01) is the component on a 
pure state of its spin-flip: 

C n (|\|r» = | (V|U|V>|/<V|V>, where 1%) = (=ia^W) = (-/o- v ) "M (2) 

The concurrence of an n qubit state with n odd vanishes identically. This monotone is noteworthy for two reasons. 
First, there is an explicit, computable closed-form expression for the minimum C„(p) which is again defined in 
terms of the eigenvalues of pp = p{a y )® n p(a y )® n (U00; BB04a). Second, in the context of concurrence dynamics 
we may study entanglement dynamics. This paper concerns itself with the latter topic, and we henceforth consider 
only pure states and unitary maps. 

The primary mathematical tool used in this paper is the Concurrence Canonical Decomposition (CCD.) This 
is discussed in detail in §11. Briefly, it is a way to decompose a unitary on n qubits into a factor that changes 
concurrence and factors that do not. Let v : ${ n — > tt n be a unitary evolution. Consider the CCD v = k\aki 
(BB04b). Now k\ and hi are symmetries of the concurrence, reducing concurrence dynamics to the second factor. 
This a factor applies relative phases to a basis of GHZ-like states. Such phases are not unique due to choices of 
diagonalization while computing the CCD, but the spectrum spec(a 2 ) is uniquely determined by v. Moreover, the 
two-qubit test for maximal entanglement capacity (ZVSW03) generalizes to n qubit concurrence capacities if n is 
even: 

Let v = k\ak2 be a CCD of v. Consider spec(a 2 ) as a subset of the unit circle. Then for n = 2p, 
there is a |\|/) e 9i n with C„(|\|/)) = and C„(v|\|/)) = 1 if and only if is within the complex hull of 
spec(a 2 ). (ZVSW03; BB04b) 

Also, for even n there is an explicit numerical algorithm for computing the CCD and hence spec(a 2 ) (BB04b). 

This work presents three new results. The first is an extension of concurrence capacities to the case n odd. For 
n even, the concurrence symmetry group K to which k\ , £2 belong is up to a similarity transform an orthogonal 
group. For n odd, K is not orthogonal but symplectic, a has repeat eigenvalues, and C2 P -i(|v|/)) = for all |\|/). 
Nonetheless, we define a two-argument lift of the usual concurrence, say C(|<|>}, |\|/)). (See Equation 7.) Suppose 
we define the amount of concurrence an odd-qubit unitary v creates to be 

k(v) =max{c(v|*),v|>|f));c(|<|>),|v))=0} (3) 
This generalized capacity has the following properties: 

• For n even, the one-argument concurrence capacity and the two-argument capacity of v coincide. 
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• For n odd, often K„(v) 7^ for the pairwise capacity despite C n (|v|/)) = 0. Further, K„(v) = 1 if and only if 
lies within the convex hull of spec(a 2 ) for any CCD by v = k\ak.2- 

• Concurrence capacity monotonicity : Using double argument capacities, the capacity of v <g> h is always at 
least that of v. 

Hence there exists a theory of odd-qubit concurrence dynamics, even though concurrence vanishes identically (on 
the diagonal) in odd qubits. 

Second, we present an explicit numerical algorithm for computing the odd-qubit CCD. Various matrix log- 
arithms must be computed, after which one invokes work in the numerical analysis literature (DGKW84) to 
diagonalize a time reversal symmetric Hamiltonian using symplectic matrices. 

We close with the third observation, which we will refer to as Kramers ' nondegeneracy : 

On the H-quantum bit state space, suppose that a O-time reversal symmetric Hamiltonian H has a 
nondegenerate eigenstate \X). Then (i) « is even and (ii) C„(\X)) = 1. In particular, \X) is entangled, 
i.e. 

The proof follows from viewing 15 as a time reversal symmetry operator in Wigner's axiomatization, a point of 
view which also simplifies the derivation of the CCD. Kramers' nondegeneracy leads one to wonder whether 
useful entangled states may be produced by cooling the system of a 15-time reversal symmetric Hamiltonian. We 
consider the perturbative stability of this entanglement while breaking the time reversal symmetry here, while the 
thermal stability of the Kramers' nondegeneracy for the Ising model is considered elsewhere (BB04a). 



II. BACKGROUND AND PRIOR WORK 

Since our key tool is a generalized Canonical Decomposition (BB04b), we review the Canonical Decomposition 
literature. The two-qubit Canonical Decomposition (CD) states that any two-quantum bit unitary evolution v, i.e. 
any 4x4 unitary matrix v, may be written: 

v = e !tp (wi ® U2)a(ii3 ® W4) (4) 

Here u\, 112,113,114 are one-qubit (2 x 2) unitary matrices, which may be chosen to have determinant one. The 
unitary a is diagonal in the Bell basis and may be thought of as applying relative phases to this basis. However, it 
is better computationally to think of a as phasing the magic basis (BDSW96; LKHC01) instead: 

|m0) = (|00) + |11))/V2, |ml) = (|01> - |10»/V2 

|m2) = (»|00>-»|11»/V5, |m3) = (i|01)+i|10))/>/5 W 

Let E be defined by E\j) = |mj), and let SU (2") denote the Lie group of determinant one 2" x 2" unitary matrices, 
50(2") denotes determinant one orthogonal matrices, and D(2") denotes the diagonal 2" x 2" unitary matrices. 
A diagonalization argument shows SU (4) = 50(4) D(A) 50(4). Moreover, the magic basis has the property that 
E'SU(2)<g)SU(2)E = 50(4), i.e. determinant one tensors have real matrix coefficients in the basis. Thus the 
canonical decomposition may be computed by transforming the diagonalization through E: 

517(4) = [£50(4)£ t ][£O(4)£ t ][£50(4)£ t ] = 51/ (2) ®SU(2) {ED{A)E 1 ) 517(2) <g) 51/(2) (6) 

We next provide a brief account and references for the best known applications and generalizations of the CD. 

Makhlin ( 102) anticipates the Canonical Decomposition by directly computing that the double cosets [SU(2) ® 
5t/(2)]\5t/(4)/[5t/(2) (£>SU(2)] are parametrized by three real parameters, the number of parameters in a given 
det(a) = 1. The CD appears explicitly in Kraus and Cirac (KC). In an important paper, Khaneja, Brockett, and 
Glaser point out that one may view the CD as an example of the G = KAK decomposition theorem for G = SU (4), 
K = SU (2) <E> 51/(2), and A = A the commutative Lie group that phases the magic (or Bell) basis (KBG01). They 
also consider the matrix factorization from the point of view of control theory in order to compute minimum times 
for applying a given two-qubit unitary evolution. Zhang, Vala, Sastry, and Whaley made use of this observation 
to describe which 4x4 unitaries v are equivalent up to tensors of one-qubit rotations. The factor a £ A is not 
unique but depends on choices of diagonalization, and these are described geometrically using Weyl chambers. 
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Specifically, the Weyl group orbit of any a produces all possible a, and each orbit intersects the Weyl chamber 
once. For G = SU (4), the Weyl chamber is a tetrahedron (ZVSW03). Finally, the terms Canonical Decomposition 
and magic basis are by now standard, and there are published surveys (e.g. (CHN03, §11. B).) The timing arguments 
of Khaneja et. al. (KBG01) have also recently been experimentally verified in liquid-state NMR (NVKT04). 

There are many applications of the two-qubit CD. In addition to timing as above, they include (i) studying the 
entanglement capacity of two-qubit operations (ZVSW03), (ii) building efficient (small) quantum circuits in two 
qubits (BM03; VW03; VD03; SBM03), and (iii) classifying which two-qubit computations require fewer than 
average multiqubit interactions (VD03; SBM03). 

Besides the CCD (BB04b), there is another n-qubit generalization of the Canonical Decomposition due to 
Khaneja and Glaser (KG01). It is also defined in terms of a G = KAK decomposition. Label N = 2" for the 
remainder. The type of a G = KAK decomposition follows from a classification theorem of Cartan involutions 
and determines the groups K and A up to Lie isomorphism. (The classification appears in Helgason (HelOl, 
pg.518). See ibid, for details.) Given G = SU(N), the three possible types demand K = SO(N) (type AI,) 
K = Sp(N/2) a symplectic group (type All,) or K = S[U(p) © U(q)] for p + q = N a block unitary (type AIII). In 
the All case, the structure of the A group also demands any a £ A has even-degenerate eigenvalues. The two-qubit 
Canonical Decomposition is type AI, and indeed the similarity transform by E shows SU (2) ® SU (2) = 50(4). 
The CCD alternates AI and All as n is even or odd. The KGD of Khaneja and Glaser technically contains two 
G = KAK decompositions, the first of which is type AIII for n > 2. In fact, the KGD is similar to the Cosine Sine 
Decomposition (CSD) of numerical linear algebra (B04) and so may be computed numerically. Physically, the 
K = S[U(N/2) ®U(N/2)] group of the KGD may be viewed as those unitaries commuting with measurements in 
the z-basis of the least significant qubit, i.e. commuting with I N i 2 ® O z . 

We next recall notation from quantum computing. The one-qubit state space is 9i\ = C{|0)} © C{|1)}. For n 
quantum bits, !H„ = = H\ ® • • • <g> !H\. (See (NCOO).) A local state |\|/) is any state which may be written 

as (g>" =1 |\|//) for € while an entangled state is any state which is not local. Notations such as e.g. |7) 
refer not to the state of a qudit but rather to a multiqubit state, e.g. |111) = |1)®|1)®|1). The n-concurrence of 
Equation 2 is an entanglement monotone (BB04b). Besides the well-known two-qubit concurrence (H97), even 
qubit concurrences n-qubits (WC01) (SC03, Eq.62) (BB04b) have also been studied. Since the single-argument 
concurrence vanishes for n odd, we introduce a two-argument generalization. 

For 15 per Equation 2, the concurrence bilinear form (BB04b) is the map C„ '■ M,, x Hn ~~ * C given by 

C„(|^>,h/>)='P%> (7) 

The complex conjugate forces the two-argument function to be complex bilinear rather than complex bi-antilinear, 
and the concurrence monotone is the norm of the form on the diagonal: C„(|(]))) = |C«(|0), |0))|. The bilinear form 
C n is symmetric for n even and antisymmetric for n odd, which causes vanishing of the monotone but not the form 
in the odd-qubit case. 

The CD is an example of the G = KAK decomposition theorem (HelOl, fhm8.6,§VIL8) for G = SU(N). This 
theorem produces a decomposition of a reductive Lie group G for any 8, a as follows: 

• The map 8 : q — > Q for g = Lie(G) is a Cartan involution (HelOl, §X.6.3,pg.518). By definition 1 , (i) 8 2 = l g 
and (ii) 8[X,T] = [QX, QY] for all X,Y £ g. As is standard, we write g = p 1 for the decomposition of g 
into the —1 and +1 eigenspace of 8. 

• Given 8, a C p is a commutative subalgebra which is maximal commutative in p. 

Note that £ is closed under the Lie bracket, while this is trivially true for a. Thus the exponential of each is a 
group. Label K = exp t, A = exp a, where for linear G C GL(n,C) the exponential may be interpreted as a matrix 
exponential. The theorem then asserts that G = KAK = {k\ak2 ; &i,&2 £ K,a g A}. 

The CD is seen to be an example as follows, cf. (KBG01). Take 8 : su(4) -> su(4) by Q(X) = 
(-/o>')^Ho>r 2 and a = span R { i|0><0| - i|l)<l| -i\2)(2\ + i\3)(3\, i\0)(3\ +i\3)(0\, i\l)(2\ + i\2)(l\ }. 
Extending these choices to n qubits produces the CCD: 



Some authors only use the term Cartan involution in the case that is a noncompact Lie algebra. In their terminology, this definition of 
a Cartan involution on the Lie algebra of a compact group, e.g. su(N), is the image of a Cartan involution of a noncompact Lie algebra 
through symmetric duality (g = t© p) < > (0 dual = t®ip). 
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Definition ELI [CCD,(BB04b)] Define 9 : su(N) -> su(N) by 0(X) = [{-ia y ) m ] f X(-ia y ) m . Then { denotes 
the +l-eigenspace of while p denotes the — 1-eigenspace. Finally, in case n is even we define 

a = S p WM ({i\j}{j\ + i\ N -j-l}(N-j-l\-i\j + l}(j + l\-i\N-j-2)(N-j-2\; 0<;<2"- 1 -2} 
U { (iV-;-l|+ j - 1) (/I ; < j < 2»~ l - 1 },) 

(8) 

with A = exp o. In case n odd, we drop the second set: 

a = span R ({ i\N - j - l)(N - j - 1 1 - i\j + 1) (j + 1 1 - i\N- j - 2) (N - j - 2| ; < ; < 2"' 1 - 2 }) 

(9) 

The Concurrence Canonical Decomposition {CCD) in n-qubits is the resulting matrix decomposition SU(N) = 
KAK. Note that n may be even or odd. 

In an earlier work (BB04b), computations in Dirac (bra-ket) notation show that Q(X) is a Cartan involution and 
a is maximal-commutative in p. The G = KAK theorem (HelOl, thm8.6,§VII.8) then shows that the CCD exists. 
Further, the CCD may be computed numerically in the even qubit case (BB04b). 

The CCD is a useful tool for studying concurrence capacities since K = exp(t) consists of symmetries of the 
concurrence form of Equation 7, where I is given per Definition II. 1 (BB04b). 

(v e K) <=> [ Cn(v\$),v\V)) = C«(|<i>>,|V» for all e 9{ n } (10) 

In particular, the above may be used to verify that SU(2) m C K as a subgroup of large codimension. One 
explanation for the fact that K alternates between orthogonal and symplectic groups is to note that the form C n is 
symmetric or antisymmetric as n is even or odd (BB04b). Another outlook, illustrated in §V, is that the spin-flip 
13 is a bosonic or fermionic time reversal symmetry operator as n is even or odd, i.e. 15~' = (— l)"U. 

III. ODD-QUBIT CONCURRENCE CAPACITIES 

The main results of this section are summarized in Theorem III. 1 1 . Each is proven in turn. 

A. Double-argument capacities generalize single-argument capacities 

To begin, we introduce a pairwise concurrence capacity K„(v) and denote earlier concurrence capacities 
(BB04b) with a tilde. 

f ic„(v) = max{C„(v| ¥ )) ; - 1,C„(| ¥ )) = } 

\k„(v) = mnx{\Cn(?W)Mv))\: <<Hf = <vlv} = i>c«(H>>,|y>) = o} {LL> 

Due to Equation 10, any CCD of a unitary v = ^10^2 implies ic„(v) = K n (a) (BB04b) and K„(v) = K„(a). 
Proposition III.l Suppose n = 2p is an even number of qubits. Then K n (v) = K„(v). 
The proof requires certain results from the literature (BB04b; ZVSW03). 

• There is an n = 2p qubit entangler £b so that for any k € K, EokE^ is a real unitary matrix, i.e. orthogonal. 
The columns of Eo resemble |GHZ) states. 

• For this Eq, any CCD v = k\ak2 moreover has d = E T aE for d = YlJ=o (j\ diagonal. As d is unitary 
diagonal, each dj is on the unit circle within C. 

• The concurrence spectrum becomes A, c (v) = {dj}^^. Then K2 n (v) = 1 if and only if £ C lies within the 
convex hull of X c (v), a subset of the unit circle (BB04b, Lem.III.2). 

• A corollary (BB04b, Scho.2.18) of the symmetry group theorem shows that £0 also translates between 
C n (-, -) and a simpler bilinear form: C„(E zi ,Eqzj) = z\z,i- 
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Example III.2 We use the CD to compute a two-qubit concurrence capacity. Consider a family of controlled- 
phase gates, e.g. v(f) =e-'' f |0)(0| +e- ! ''|l)(l| + e- ,Y |2)(2| + e 3,Y |3}(3| with det[v(f)] = 1. ApossibleCD is: 

v(r) = (e-'''° : <g> I 2 ) e itaZmZ (h ® e"^) (12) 

The central factor is a valid choice for a in v(t) = kiak2, since e !fa ®° is also diagonal in the magic basis. Thus 
X c [v(t)] = spec(e 2i ' mW ) = {e 2,r ,e 2 ''',e- 2, '',e- 2 '' f }. Only for t G |Z do we have within the convex hull of X c [v(t)], 
and the convex hull theorem asserts iC2[v(7t/4)] = 1. Indeed, up to phase v(tc/4) = |0)(0| + |1)(1| + |2)(2| — 

|3)(3|. Moreover, if H = ^ j j ^ is the Hadamard gate (NCOO), a standard identity converts v(k/A) into 

the quantum controlled-not: 

CN0T= |00)(00| + |01)(01| + |10)(11| + |11)(10|= {h®H)v(%/A){I 2 ®H) (13) 

Thus v(;t/4) carries an unentangled state to a maximally entangled state, since CN0T(H® I 2 ) |00) = CNOT-^ (|00) + 
1 10)) = ^(|00) + |11)). More intricate examples intwo-qubits (ZVSW03; SBM03) and an even number of qubits 
(BB04c) are available in the literature. 

Lemma III.3 Suppose the number of qubits is even. Let z\ = T!J=o a j\j)> Z2 = T!J=o bj\j), and zi = T!J=o c j\j) 
throughout, and let X c {v) = {Xj}^^. Then we have the following. 

K„ (v) = max { | L^-l c 2 Xj | . ,t Z3 = j > ^ = Q} 

K„ (v) = max { | ijjo 1 A,; | ; z\zi = z\zi = 1 , z\zi = 0} 

Proof of Lemma III.3: The first equation appears in (BB04b); cf. (ZVSW03). For the second, take vectors 
z\ , Z2 and label x = Eqzi , y = EqZ2 ■ Then 

[C n (x,y) = 0] [c„(£ozi ,£oz 2 ) = 0] [z\z 2 = 0] (15) 

Moreover, without loss of generality by choice of z\, Z2, and symmetry we may suppose v = E^dE^ for d 2 = 
l!jZo h\J) Ul Then C n {EodE^x,E dE^y) = C„{Eodzi,Eodz 2 ) = {z\d T )dz 2 = E^, 1 a J b J X P D 
Proof of Proposition III.l: Let aj, bj be chosen so as to maximize the expression for K„(v) per Lemma III. 3, 
i.e. K„(v) = | L^=o a jbjXj\ - Now choose complex numbers cj so that c 2 = a jbj, and put zi = Lj=o c j\j)- We note 
that z\ Z3 = 0. Moreover, z^zj < 1, for 



N-l N-l N-l N-l i i 

E \cj\ 2 = I |c?|= £ \ajbj\ < £ -|a/ + -|fc/ = l 



(16) 



j=Q j=o j=o j=o 
Label t 2 = Z3Z3, noting t 2 < 1. Then (f _1 Z3)^(f _1 Z3) = 1, so by definition of K2 P (v) we have 

iV-l 7V-1 
K„(v) > K„(V) > I £ r 2 C 2 ?l 7 -| = r 2 | £ flyiyXyl = r- 2 K„(v) (17) 

y'=o j=o 

Thus r = 1 and hence K„(v) = K„(v). □ 



B. Monotonicity 

We next demonstrate concurrence capacity monotonicity, i.e. that j i— > K n+ j(v<g> if 1 ) is monotonic. 

Proposition III.4 Let n be either even or odd, v € SU(N) an n-qubit computation, and let I2 denote the trivial 
one-qubit computation. Then VL n+ \(y®l2) > K„(v). 
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Proof: Choose \\\t) such that k„(v) = C„(v|$),v|\|/}) while C„(|<|>), |\|/)) = 0. Then |(j)) ® |0) and |\|/) (g) |1) are 
a null-concurrent pair of (n + l)-qubit states: 

C„+i(|4)>®|0),| ¥ >®|1>) = (^|®(o|)Ha>r n+1 (l¥)®|i» = [c„(|^>,| ¥ »](<o|Ha- v )|i)) (18) 

Now (0|(— /a J )|l)) = 1, so the above expression is [0](1) =0. A similar argument demonstrates that 

C„ + i[(v®/2)(|^)®|0)),(v®/ 2 )(| ¥ )®|1))] = [C„(v|4>),v|\|/))] [Ci(|0),|l))] (19) 

The second term of the product is one, while the first is K„(v). Since we have exibited a pair for which v<E>h raises 
the pairwise concurrence by at least K„ ( v) . Since K n+ 1 (v Cg> h) is the maximum over all null-concurrent pairs, while 
\§) ig) |0), |\|/) <g> 1 1 ) is such, we see K„ +i (v® 7 2 ) > K„(v). □ 



C. Parity-independent concurrence spectra 

We extend the maximal concurrence capacity condition of Zhang et. al. and Bullock, Brennen (BB04b) to 
odd-qubit systems. The first step is a definition valid in either parity. 

Definition III.5 Let v £ SU(N), N = 2". For n of either parity, the concurrence spectrum X c (v) is the set 
A- C (v) = spec ia y )® n ]^v(— ia } )®"v T ^j . Viewing v as an R-linear map, equivalently X c (v) = spec (yl5v T l5~ l ). 

We briefly show this coincides with the definition of the even-qubit concurrence spectrum of the literature 
(BB04b). The definition ibid, states that the concurrence spectrum is the spectrum of (EqvEo)(Eq\'Eo) t ■ Indeed, 
given EqEI = {-ia y ) m per the classification of £ with ESO(N)E^ = K ibid., 

spec (ElvEo){ElvEo) T = spec(£([v£o£o v r £b) = spec[(E El) f vE E^v T ] = spec[(-ia y )® n v{-ia y )®"v T ] 

(20) 

In fact, the same argument shows that X c (v) is the spectrum (E^vE)(E^vE) T for any E as above. Cf. (M02). 

The odd-qubit case requires different similarity matrices, say F (BB04b), which translate K not into an or- 
thogonal group but rather a symplectic group per Equation 22. For the concurrence form C n (— , — ) for n-odd is 
antisymmetric, and symplectic rather than orthogonal groups are the appropriate symmetries of antisymmetric 
bilinear forms (i.e. two-forms.) For a standard similarity matrix, we take 



= I^o~' \J) 01 + \N~j- 1) 01 + (N/2 + j\ -\N-j- l)(N/2 +j\), where 
{l.,-}^- 1 C {±1} by (-^r =lJio"\(PV-J- 1)01 - \J)(N-j-l\) 

Also, label throughout Jn = {—io y ) ®In/2- Before showing that Fq translates K into the standard symplectic 
group, we show that Fo carries c(— , — ) to the standard two-form J4 ( — , — ). 

Lemma III.6 For A (|(j)), for)) = W\Mv)> C„{Fo\ty},Fo\y})=A(\ty},\y})forall\§},\\\f) e 

Proof: C„(F \§),F \\\r)) = W\ F o (-io y )® n F \\\r). Now F J N Fj = (-io y )® n (BB04b, PropII.14), whence 
FT{-ia y )® n F = J N . □ 

Now Sp(N/2) is that copy of the symplectic group which embeds within SU (N) as the symmetries of A (— , — ), 
i.e. satisfying A (v|0),v|\|/)) = A. (|<|)}, |\|/)) for all G 0i n . In block form: 

S P (N/2) = { V eSU(N);v T J NV = J N }=U* *)eSU(N); ^C^y^Mc,^ is symmetric, 1 (22) 
As Eq SO(N) El = K 2p , so too Fq Sp(N/2) Fj = K 2p -i. 



8 



We next associate X c (v) to spec(a 2 ) for v = k\ak 2 in the odd-qubit case. Suppose we label D to be the following 
diagonal subalgebra of SU (N): 

{N/2-1 N/2-1 "I 

E dj(\j)(j\ + \N/2 + j)(N/2 + j\); n dj = ±l (23) 
j=o j=o J 

Now there is a standard SU(N) = KAK decomposition which follows from 8aii(///) = JN{—iH T )jl/ (HelOl, 
§X.2,pg.452) and a = logD as above. Given a v G SU(2 2p ~ l ), it writes v = coi^G^, with coy G Sp(N/2), j = 1,2 
and d eD. 

Suppose given v G SU(2 2p ~ l ), we then write FqvFq = 03ido) 2 , with coy G Sp(N/2), j = 1,2 and d £ D. The 
odd-qubit CCD again follows by a similarity transform: v = (Fq(£)iFq ){Fq<1Fq )(Fo(02Fq ) with a = FqcLFq G A, 
£; = Fq&jFq G /T, y = 1,2 is a CCD. Note that a is diagonal on the GHZ-like basis states {fob')}- 

Lemma III.7 Lef n = 2p — 1. Then for v = (Fo<»iF r ) (FodF r ) (Fo<»2F(f ) fne CCD as above w;f/i if = 

L^io 1 dj(\j) (y | + |7V/2 + j) (N/2 + y'|) diagonal and determinant one, we have X c (v) = {d 2 }^/^ 1 U {d 2 }^/^ 1 
(counted with multiplicity.) 

Proof: Given A,B, invertible, spec(AB) = spec(BA). Also, Equation 10 is equivalent to a matrix equation 
k r (-iti*)® n k = (-io*)®" for all keK. Recall FqJ n FJ = (-/a- v )®". Then 

X c (v) = spec^[(-/o- v )®"] t v(-/o- v )®"v :r j = spec^[(-/a- v )®"] t fcifl/t 2 (-/o>') 8 "/t 2 r fl 7 '/t[ 
= svec(k\[{~ioy) m ]k x ak 2 [(-i&Y m ] T k T 2 a T ) = ^(^[(-^)^]Jfcia[j£(-^)^Jfc 2 ] r a r 



(24) 

spec^[(-^)®"H(-i^®"] r a r J = spec^- [FoJ N F^]F dFj [FoJ N Fj]F d T Fj j 
= spec( —FoJ^dJ^d 7 Fq ) = spec( —J^dJ^d ) = spec(c/ 2 ) 

The last equality makes use of d G £> repeat diagonal. □ 
D. A convex hull argument in odd qubits 

Definition III.8 Suppose n = 2p— 1. The reduced concurrence spectrum X c (v) of v £ S£/ (iV) is the set {A,y}^Q 1 
for v = fei (FodFQ)k2 a canonical decomposition of v and c/ = LjZo 1 v^7(U) 01 + W/2 + j) (N/2 + y'|). The 
convex /iw/Z CH[X e (v)] of X c (v) is the set convex linear combinations of the points of X c (v), i.e. 

{N/2-1 N/2-1 ~| 

E 'A/ ; < < 1 , E *) = hh* Mv) f ( 25 ) 
7=0 >=o J 

Proposition III.9 Suppose n~2p—\ is an odd number of qubits. Throughout, label zi = Ly=o a ;'l./)> ~2 = 
Y!J=q b j\j), and X c (y) = {Xj}^/^ Then the following hold. 



K„(v) = max 



ifjll 1 h ( a N/2+ jbj - a Av/2+ y) 



; z[-/iVZ2 = 0, z\zi = Z 2 Z2 = 1 



(K»(v) = l)^(0 6CHp e (v)] ). 



Proof: The first item follows from Lemma III. 6, substituting x = Fqzi, y = FqZ2- We continue to the next item. 
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For the second item, we first prove =>. If K„(v) = 1, then we may choose zi, z 2 so that 

1 = \l!Jll~ l ^j( a N/2+jbj-ajb N/2+ j)\ < L% 2 Q ~ l \a N /2+jbj-ajb N/2+j \ 

< L^^^ + Wn/z^H 2 + \b m+j \ 2 < 1 

Here, note that the second inequality is an iterate of C\ (10) , |\|/) ) < y/ (010) (vlv), f° r a U 10) 1 I V) G #i ■ The last 
inequality in Equation 26 is the Schwarz inequality. 

Now label ay = a N / 2 +jbj — ajb N / 2 +j, for < j < N/2 — 1. Then by Equation 26, 

N/2-1 N/2-1 N/2-1 

i = l E V*;l= E IV*;I= E N ( 2? ) 

j=0 7=0 7=0 

Thus there must exist some z S C, zz = 1, so that XyOCy = z|(Xy|, and moreover L^io 1 I 01 )! = O n the other hand, 
z\JnZ2 = demands that = Ly^o ' a y = z ^j=Q ' l a /l^-/- Multiplying by z and taking the complex conjugate, 
= Ly'Zo 1 \ a j\^j which given L^fo ' \ a j\ = 1 by Equation 27 demands G CH[X, c (v)]. 

Consider now the converse case, i.e. € CH[X, c (v)]. Then there exist fy real, nonnegative so that = L^io ' tft^j- 
For < j < N/2 — 1, label complex numbers a,- = tjkj, so that we have 1 = L^io ' l a 7'l an d moreover = = 
^=o 1 0^-7 = ^J=o * a j- ^ e are re d uce d to the following question: May we choose {fly} w) 1 , {^y}y^To so mat 

AT-l AT/2-1 

a J = a N/2+jbj-ajb N/2+j and E l fl / = E N 2 = 1 ( 28 ) 

7=0 7=0 

To do this, write ay = |oty |e ,arg a J, and take ay = yjayj, a-jv/2+y = 0, bj = 0, and b N / 2+j = -e ,arg a J \f\Oj\. Then 
we see that a N / 2+ jbj — ajb N j 2+ j = ay. Moreover, 

/V/2-1 

|fly| 2 + K /2+ y| 2 = |ay|, N 2 + |^/ 2+ y| 2 = l«yl and £ |ay| = 1 (29) 

7=0 

Thus the vectors zi, Z2 per the statement of the proposition are normalized to be norm one. □ 

Hence, a similar picture emerges for concurrence capacity one for the even and odd-qubit cases with these 
definitions. The new feature, doubly-degenerate eigenvalues in X c (v) arising from the D above required for type 
All, will a posteriori be an instance of Kramers' degeneracy; see §V. 

Corollary 111.10 For n = 2p - 1, lim pMao da( { a e A ; K„(fl) = 1 }) = 1. 

The proof of the Corollary follows by considering probability density functions on the unit circle (BB04b), 
given that the number of concurrence eigenvalues grows exponentially with n. Thus most unitary evolutions 
for large « (of either parity) are maximally entangling as measured by concurrence. It would be interesting but 
technically challenging to restate this in terms of Haar measure du on SU(N). The difficulty is that the pullback 
measure from the K x A x K to SU(N) is singular, namely singular near the set where the A factor is an identity. 
For future reference, we summarize the concurrence capacity results of this section. 



10 



Theorem III.ll Let K„(v), K„(v) be the pairwise concurrence capacity and concurrence capacity respectively. 

1. The pairwise capacity and the capacity are equal in any even number of qubits. Thus, 

2. For n either even or odd, any CCD by v = k\aki satisfies K„(v = k\akj) = K„(a). 

3. For any n, we must have K„ + i (v ® I2) > K„ (v). 

4. Suppose n — 2p — 1 is odd. Then for da the Haar measure on A, 

lim Prob(x„(a) = 1 ) = lim da({a G A ; K„(a) = 1}) = 1 (31) 



IV. AN ALGORITHM COMPUTING THE ODD-QUBIT CCD 

In this section, we close a gap in the literature. Specifically, we present an algorithm for computing the CCD 
when the number of qubits is odd. We make use of an algorithm (DGKW84) by Dongarra, Gabriel, Koelling, and 
Wilkinson cited in a survey (BBM92) of diagonalization arguments. The algorithm (DGKW84), which appears 
in the numerical matrix analysis literature, improves the numerical stability and computational efficiency of the 
earlier work on time reversal by Dyson (D61). 

Recall from §111. C that it suffices to compute the standard type All KAK decomposition given by SU (N) = 
Sp(N/2) D Sp(N/2) with D the repeat diagonal subgroup of SU (N). Forgiven v G SU(2 2 ''- 1 ) for which we wish 
to compute the CCD, suppose we obtain FqvFq = G)id(£>2, with co ; - G Sp(N/2), j = 1,2 and d E D. Then v will 
have CCD v = k x ak 2 = {F (OiF^)(F dF^)(F O)zFQ). Before computing SU(N) = Sp(N/2) D Sp(N/2), we make 
one new definition. 

Definition IV.l Let H G C NxN be Hermitian. Recall J N = (-ia y ) <E> I N / 2 . We say that the Hamiltonian H is 

J^-skew symmetric iff HJn — JnH t = 0. 

Remark IV.2 In (DGKW84), the above is the definition of "H has a time reversal symmetry." Indeed, time 
reversal symmetry follows for the operator © = J^x, (x complex conjugation,) per the upcoming Definition V. 1. 
Moreover, for the standard type All Cartan involution (HelOl, pg.452) 9aii(X) = J N XjJj, let su(A^) = Paii© ^aii 
for the corresponding Cartan decomposition into —1 and +1 eigenspaces. Then H is Jn skew-symmetric iff. 
iH G Paii- Indeed on su(N), X = —X T . Hence — iH = J^WJ^ = -JnIH 7 ^ if and only if HJ N = J N H T . 



A. Algorithm for the standard All KAK decomposition, SU(N) = Sp(N/2) D Sp(N/2): 

The outline below for computing the standard SU (N) = KAK decomposition of type All (see §111. C) is similar 
to the AI case used in (BB04b) to compute the even-qubit CCD. The added difficulties are (i) a more complicated 
formula for p 2 and (ii) a more delicate diagonalization argument for p 2 once computed. In fact, the latter requires 
the symplectic diagonalization argument referenced above. 

Lemma IV.3 Suppose v G SU(N) with v = pk for p = exp(iH) with H a Jn skew-symmetric Hamiltonian and 
k G Sp(N/2). Then p 2 = -vJ N v T J N . 

Proof: We have H T = —JnHJn, given = jjj = —J^. Thus for any t G E, [exp(iHt)] T = jjjexp(iHt)JN = 
—JNSxp(iHt)JN This holds in particular for p. Now put w = v\ so that w = kp for k = k\ p = p^ . Thus 
p T = jljpJfj. Moreover, k G Sp(N/2) demands k T J^k = J^, as Sp(N/2) is a group. Thus —J^w T J^w = p 2 . 
Taking the adjoint of each side produces the result. □ 

With this lemma, we now present the algorithm for computing the standard type All decomposition. 



11 



1. Suppose v = pk per Lemma IV.3. Compute p 2 = —vJnv t Jn. 

2. We may write p = exp(iH) for some Jn skew-symmetric Hamiltonian H. Compute a logarithm of p 2 = 
exp(2iH). The diagonalizing matrix implicit in computing the matrix log need not be symplectic, and 
generic logarithms will take the form 2iH for some (2)H which is Jn skew-symmetric. 

3. Compute a symplectic matrix tOi € Sp(N/2) so that 1H2 — <X)\(iH)(Oi is repeat diagonal, per §IV.B. 

4. Label p = CDiexp(;7/2)coj and d = exp(;7/2)- Compute 0)3 = p^v. Then CO3 G Sp(N/2). 

5. Put CO2 = (o[m G Sp(N/2). Note that coififcoj = p. Thus the type All decomposition is v = [001] [d] [cojab] = 

(8\d(£>2 

This concludes the overview of computing SU (N) = Sp(N/2) D Sp(N/2). The next subsection details Step 3. 



B. Symplectic diagonalization 

In this section we address the problem of finding the eigendecomposition of a matrix H which is Jfj skew- 
symmetric. Generically, these techniques work on any square matrix with an even number of rows and columns, 
and there are no simplifications when the size is a power of two. Thus we describe the generic case where 

Jit = ^ ^ q J and H = is also Jq£ skew symmetric. 

( A B \ x T 

Explicitly, J^-skew symmetric means H = I — -j J , where A = A' and B = —B are t x I matrices. We 

will construct a unitary skew-symmetric Hamiltonian matrix co of the form co = ^ ^ ^ ^ , so that the columns 

of co are the (right) eigenvectors of H. Each eigenvalue A* for k = 1, . . . ,£ of H is real and of multiplicity 2. In 
particular, both the kih and the (£ + k)^ 1 columns of CO are eigenvectors of H corresponding to A,^. Also, given the 
block form, co E Sp(N/2) up to global phase. 

The algorithm of Dongarra et. al. (DGKW84) proceeds in two major steps. First we reduce H to block diagonal 
form using a similarity transformation, and then we use the QR algorithm to find the eigenvalues of the blocks. 
We consider each of these phases in turn. 

First, we construct a skew-symmetric Hamiltonian unitary matrix Q of the form Q = 

^ ) so that QHQ? = f « T* ) where T is real, symmetric, and tridiagonal. We initialize 
Q2 Qi J V u 1 ) 

Q to be the 21 x 2£ identity matrix. In order to preserve the structure, we construct Q as the product of two simple 
types of matrices: 



ci I r — b i y 

and if we let r L = \a\ z + \b\ 2 , then a matrix of the form ( -r . I is unitary. In addition, 



The product of 2 x 2 skew-symmetric Hamiltonian matrices is also skew-symmetric Hamiltonian, 

a/r ~b/r 
b/r a/r 

air -b/r \ f a b\ / r 0\ 

r , , r — = n so the unitary matrix can be used to introduce zeros. Choose / be- 

b/r a/r J \ -b a J \0 r J 1 J 

tween 1 and i and construct a matrix R as the 2£ x 21 identity matrix except that entries Re+j.t+j = Rj / = a/r 

and Rj.i+j = —Rt+j £+j = —b/r. Then the product RH is equal to H except that the entries in rows j and 

C + j become 



(RH)j. k (RH) j/+ k \ = (l/ r - b / r \ ( A J^ hi 
{RH) i+jik (RH) t+j>l+k ) \ b/r a/r ) \ -B } , k A uk 



(32) 



k = !,...,£. Since this product is skew-symmetric Hamiltonian, so is RH, and it can be shown in a similar 
way that (RH)R 1 ' is skew-symmetric Hamiltonian. Thus we can use R as a similarity transformation that 
preserves the structure. 
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• Let S be a real orthogonal matrix of dimension £ X t. Then 

(S 0\ / A B\(tf W SAS< SBS<\ 

[o S ){-B A ){o S+ )-\ -SBS^ SAS^ ) (33) 
is skew-symmetric Hamiltonian. 

Using these matrices, our construction takes £ — 1 steps. We describe the first step in detail. 

The first step places zeros in the first column of the matrix in rows 3 through 21. To put a zero in position 
(£ + j, 1) (J = l,...,n), we construct an R matrix involving rows j and I + j. If r 2 - = |A ; -,i | 2 + | 2 , then this 
matrix Rj is the identity matrix except that entries R( + jj + j = Rjj =Aj^/rj and Rj.e+j = —Re + jj + j = —Bj^/rj. 
We replace H by (RH)R 1 ' and update Q by premultiplying by Rj, repeating this for j = 1 , . . . , I. 

We complete the first step by putting zeros in rows 3 through £ of column 1. Note that these elements are 
now real, since elements 2 through £ are just the values ry Thus we can construct a real orthogonal reflection 
(Householder) matrix of the form S — I — 2ss T where § = [0,r2 + \\r\\ ,n, . . . ,r„] T and s = s/\\s\\. A similarity 

transformation of H by ( ^ ^ J produces the required zeros, and Q is updated by premultiplying by this matrix. 

Steps 2 through £ — 1 are similar; in step k we first put zeros in the B portion of column k using R matrices and 
then zero elements k + 2 through £ of the A portion using a reflection matrix. The final result is that the transformed 
H has a real tridiagonal matrix T in place of A and A and zeros elsewhere. 

The QR algorithm is considered to be the algorithm of choice for determining all of the eigenvalues and eigen- 
vectors of a real symmetric tridiagonal matrix. We use the algorithm to form X, the matrix of eigenvectors of T. 
Implementation of the algorithm requires care, and high quality implementations are available, for example, in 
LAPack (LAPACK). Other codes are available at http : / /www . net lib . org. 

We construct the eigenvector matrices U and V as U = Q[X and V = Q-^X. Note that most implementations 
of the QR algorithm do not guarantee that the eigenvalues are ordered, so a final sort of the eigenvalues and the 
columns of U and V should be done at the end if desired. 



V. TIME REVERSAL, THE CCD, AND KRAMERS' NONDEGENERACY 

The section presents three topics, all following from an interpretation of O from Equation 2 as a time reversal 
symmetry operator. First, the Cartan involution defining the CCD may be rewritten entirely in terms of the spin- 
flip, and the eigenspaces of Q(iH) are associated to time symmetric and antisymmetric Hamiltonians H in a natural 
way. Second, a well-known procedure exists to convert any G = KAK decomposition into a polar decomposition, 
and the polar decomposition associated to the CCD writes a unitary v € SU (N) as a product of two factors, one 
evolution by a time symmetric Hamiltonian and one evolution by a time anti-symmetric Hamiltonian. Third, we 
demonstrate the entangled eigenstates of Kramers' nondegeneracy as described in the introduction and consider 
the perturbative stability of this entanglement under time reversal symmetry breaking. 



A. Spin-flips as time reversal symmetry operators 

Recall the Bloch sphere (e.g. (NCOO)), which provides a picture of the data space of one qubit. As a remark, 
the Bloch sphere may be thought of as a parameterization of the complex projective line CP 1 (e.g. (M84, §40).) 
Briefly, CP is the set of all equivalence classes of vectors in C 2 up to multiple by a nonzero complex scalar. To 
associate such a class with a Bloch vector, normalize | \|/) as above so as to write | \|/) = re" [cos 1 1 0) + e ,<p sin 1 1 1 ) ] . 
The Bloch sphere vector of \\\t), say [|\|/)],is given in spherical coordinates by (1,0, cp) (NCOO, pg.15). Recall also 
that the north pole is [|0)] and [|1)] is the south pole. 

Now let b G (F2)' 1 be an ra-bit string. The typical procedure when quantizing a classical computation is to 
extend the classical outputs linearly without phases. Thus, a reasonable interpretation of quantum bit-flip would 
be (cr*)®". This is the common interpretation, but note that in one qubit a x is not reflection on the Bloch sphere 
and indeed has a fixed state, (1/V2)(|0) + |1}). Rather, the odd reflection of a single qubit under the Bloch 
parameterization of CP 1 is the spin-flip \\\f) h-> (— z'<j v )|\|/) = (— ;'o- v )|\(/). 
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The appropriate physical interpretation of the spin-flip is as a time reversal symmetry operator (W59, Ch.26) 
(G66, pp.3 14-322) (K64; S85). Wigner defined a generic time reversal symmetry operator © as any M-linear 
involutive map of the quantum Hilbert space which is antiunitary, i.e. complex anti-linear (@(a|\|/i) + PIV2}) = 
a©|\|/i) + p®|\(/2)) and orthogonal in the induced real inner-product on K 2 '' = C p . Generic time reversal symmetry 
operators are usually denoted by a capital ©; we ask the reader's forebearance in distinguishing this from the 
lower-case describing a Cartan involution. 

Such a time reversal symmetry operator ® maps the state of a system to its motion-reversed state, so that 
momentum eigenstates transform as ©|p) = |— p). In particular, if our qubit is a spin j particle, e.g. with |0) = 
I |) and |1) = I I), then \\\t) G 5{\ rotates counterclockwise about the Bloch sphere vector of [|\|/)]. Thus 15 
per Equation 2 is the natural quantum angular momentum reversal in n-qubits. Indeed, the total spin angular 
momenum, S = Y!j=\Gj> is inverted under time reversal: USU^ 1 = —S. Spin-flip operators may be defined for 
(f-level systems (qudits) but may not both preserve pure states and commute with local unitaries (RBCHM01). 

We note in passing that the spin-flip picture also allows one to quickly rederive one of the monotone properties. 
Namely, antipodal points in the Bloch sphere parameterization of the complex projective line CP 1 correspond to 
Hermitian-orthogonal states of H\. Hence, C„(|\|/)) = | (\|/|15|\|/) | — if |\|f) =g)y =1 |V;) (the monotone property,) 
since in this event has a factor (Yyi^lV;) = 0. More generally C„(|\|/)) = whenever \\f) = \\f\) ® ^sf-i) 

for G ^,,-1 and [|\|/2)] a point on the Bloch sphere. However, the latter is not an equivalence for n even. 
Consider |W 4 ) = (l/2)(|000lj ) + |0010) + |0100) + |1000». 



B. Time reversal and the CCD Cartan Involution 

We next show that physically, the eigenspaces of the Cartan involution producing the CCD correspond to 15- 
time symmetric and IS-time antisymmetric Hamiltonians. They are then explicitly described in the Pauli-tensor 
basis of su(N) in much more compact form than in Dirac notation (BB04b). 

Definition V.l Consider H a Hamiltonian on a finite dimensional Hilbert space tt\ i.e. H is selfadjoint within 
Endc(^) C Endjj(.?f ). Then H is time reversal symmetric with respect to iff H = &H© 1 as elements of 
Endij(^). A Hamiltonian is time reversal anti-symmetric with respect to © iff // — — ©H® \ 

Proposition V.2 LetQ(X) per Definition ILL Label su(N) ~ p@i as the —1 and+l-eigenspacesofQ. LetUbe 
the spin-flip. Then (i)forH a traceless Hamiltonian, so that iH G Su(N), Q(iH) = 15 (iH) I5~ l , with the right-hand 
side viewed as a composition of ^-linear maps. Also (ii) (H has time reversal symmetry with respect to 13) 
(iH G p, ) and (Hi) (H has time reversal anti-symmetry with respect to 13) (iH G t). 

Proof: Let x denote the complex conjugation operator i— ► Then 13 = (— ;a y ) 8 "x = x(— ;'a- v )®", given — ia y 
real. So IT 1 = x[(-;a- v ) g) "] t . Moreover, [(-ia y ) m ] f = (-L N ) n (-ia y ) m . Finally, z(iH)z = IE. Thus, 

U(iH)l5- 1 = (-io y f"T(iH)x[(~io y f"}^ = (-I N ) n (-ia y f n JiH)[(-ia y f n = Q(iH) (34) 

The latter two items follow at once. □ 
With the above proposition, we may describe the infinitesimal Cartan decomposition su(«) = p © i directly in 
terms of tensors of Pauli operators. Let j denote either 0, x, y, or z, with a' = h in case j = and Pauli matrices 
a*, O y , or a z as appropriate otherwise. A multiindex / = j\ ]2 ■■■ jk- • • jn denotes a string of «, and J will be said 
to be nonzero if some j k ^ 0. Finally, let /a® 7 denote z'®£ =1 (o jk ). Then su(N) = © a n non zero /R{i(J®' / } We have 
the following corollary, discovered independently by Bremner et. al. (BDNB04, Thm.5). 

Coroiiary V.3 Continue the convention of the previous paragraph, and write 

bu(N)=( M{io 8/ })©( R{i<P>}) (35) 

The above is the infinitesimal Cartan decomposition of ' Q(iH), i.e. p = ®#j = o mo d2^-{^' J }j an d % = 
©#y=i mod2^-{'^' J } I n particular, K is the Lie group of those unitaries which are exponentials of Hamiltonians 
with time reversal anti- symmetry w.r.t. 13. 
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Proof: Distinct Pauli matrices anti-commute, each has (o 7 ) 2 = h, and a y is purely imaginary while a x , O z , and 
I2 are real. Considering the tensors case by case completes the proof. □ 



C. A time reversal polar decomposition 

We next consider the polar decomposition which may be derived from the CCD. In most treatments, the polar 
decomposition of a general Cartan involution is proven and then aG = KAK theorem is derived from it. We 
next use the CCD to produce a polar decomposition for time reversal symmetry. This practical decision avoids 
rearguing the G = KAK theorem for compact groups (HelOl, thm8.6,§VII.8). 

Corollary V.4 Suppose v £ SU (N) is a phase normalized quantum computation in n qubits. Then we may write 
v = exp(iHp)exp(iH() for some Hamiltonians Hp, Hi such that Hp has time reversal symmetry and H{ has time 
reversal anti-symmetry with respect to the spin-flip 15. 

Proof: Let v = k\aki be the CCD of v £ SU(N). Then in particular v = {k\ak\)(k\k2). Since K is a group, 
k\ki is a time antisymmetric evolution by Proposition V.2. Moreover, let a = exp iH for iH £ a C p a time 
symmetric Hamiltonian. As iH £ p, we have Q(iH) = [(-ia y ) m ^ (iH)(-i&')® n = iH. Moreover, k £ K is 
a symmetry of the concurrence form (Eq. 10,) which as a matrix equation demands k 1 ' (—ia y )®"k = (—iG y ) m . 
Hence kj(-ia y f m = (-ia y ) m k\, and for k\iHk\ £ p: 

Q(hiHk\) = [{-ia y )® n ]%TiH)k\{-i(5 y )® n = ki [(-ia y )® n ] f Wj(-ia y f n kt = -h {M)k\ (36) 

Thus k\ (iH)k\ has time reversal symmetry, and the usual matrix exponential formula (valid since SU (N) is linear) 
shows k\ak\ = 

Remark V.5 Note that the vector space decomposition su(A^) = p © t makes clear any such v may be ap- 
proximated by rapid pulsing of the time symmetric and anti-symmetric factors, by applying the Trotter formula 
(e.g. (NC00, §4.7.2).) However, the decomposition above requires no such pulsing of the time-symmetric and 
time-antisymmetric Hamiltonians. 



D. Kramers' nondegeneracy 

Finally, we rederive Kramers' degeneracy in the case of 15 and note a further, 15-specific nondegeneracy prop- 
erty. Recall Kramers' degeneracy (K30; K64) proves that the eigenstates of a collection of an odd number of spin 
i electrons become doubly degenerate in the exclusive presence of a time-symmetric interaction, such as an elec- 
tric field. The degeneracy is broken with the introduction of a magnetic field. In terms of an energy Hamiltonian 
H of the system, the degeneracy corresponds to 2 or greater dimensional eigenspace for energy eigenstates. 

Lemma V.6 Suppose that £ J{ n is an eigenstate of some traceless Hamiltonian H which has time reversal 
symmetry, with eigenvalue A, £ R. Then the spin-flip 15|\|/) is also an eigenstate of eigenvalue A,. 

Proof: Since iH has time reversal symmetry, Q{iH) = -iH. Thus {-ia y ) m (iH) +JiH){-io y )< s >" = 0, and taking 
a complex conjugate produces (-ia y ) m (iH) + (iH)(—iG y )® n = 0. Now (iH)\y) = A,|\|/), so that 

(iff)U|y) = (M)(-ia y )® n '\yj = -(-iG y )® n (iH)\y) = -(-ia y f n 0^) = (37) 
This concludes the proof. □ 

Theorem V.7 (Cf. Kramers' degeneracy, (K30; K64)(S85, pg.281)) Let H be a traceless Hamiltonian on some 
number n of quantum-bits. Suppose H has time reversal symmetry with respect to 15. Let A, be a fixed eigenvalue 
of H. Then either (i) X is degenerate with even multiplicity or (ii) the normalized eigenstate \X) has C„(\X)) = 1. 
For n odd, case (i) holds: all X are degenerate with even multiplicity. 
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Proof: Let Xj be some eigenvalue of H. By Lemma V.6, both \Xj) and l5\Xj) are energy eigenstates. Should 
these two states be linearly independent, then Xj is degenerate. If any eigenvalue is non-degenerate, say X^, then 
by antiunitarity of U, we must have lS|A.jt) = e'f |A,^) for some global phase (p. Using C„ (|A,^)) = KA^UIA.^)! we 
see that this eigenstate must have concurrence one. 

Suppose in particular n — 2p — 1. Then C n (— ,— ) is antisymmetric and vanishes on the diagonal, implying 
(Xj\l5\Xj) = for all j. Consequently, \Xj) and U\Xj) are Hermitian orthogonal and may not be dependent, 
implying case (i). □ 

Thus, for the spin-flip 13 there is in addition to the Kramers' degeneracy a Kramers' nondegeneracy. As always, 
if n is odd so that the total n-qubit system is a fermion, then a time reversal symmetric Hamiltonian implies that all 
energy eigenstates are degenerate. Yet moreover in the specific case of O and n even, a nondegenerate eigenstate 
must also have maximal concurrence and hence be entangled. 

We provide some illustrative examples. First note that there are many systems endowed with time reversal 
symmetric Hamiltonians. In particular, any system with (exclusively) pairwise nearest neighbor coupling between 
qubits has iH G p, by Corollary V.3. An example of an interaction that occurs in many solid state systems is the 
quantum XYZ model: 

H X yz= £ ^(fjal+Jydjal+J^jal (38) 

<j,k> 

with J x j_ z £ R where the sum is taken over all nearest neighbor pairs and the boundaries may be fixed or periodic. 
In one dimension, these nearest neighbor coupled systems are known as spin chains. Spin chain Hamiltonians 
are of great theoretical interest, for under the appropriate parameter regime they exhibit long range classical 
correlations near a quantum phase transition (L61). We can characterize the dynamics of entanglement in spin 
chains using the concurrence capacity. With this goal in mind we observe the following useful fact: 

Proposition V.8 Let p, t be as in Corollary V.3. IfiH € p and H e R NxN , then X c {u = e - iHt ) = {er 2, V} where 
(€f parametrizes time and Xj G R are the eigenvalues ofH. 



Proof: By Definition III. 5 the concurrence spectrum of the unitary generated by iH, u = e lHt is 

X c {u) = s V ec[{-iGy)®rte- iHt (-iay)® n (e- iHt ) T ] = spec( e -'' // ' e -'' //r ') 

= spec(e- 2iHt ) = { £>~ 2 'V ; Xj 6 spec(tf) } 



(39) 



We have used (~ia y ) m UH(-ia y ) m = -iH and therefore (-Wf^Hi-io^®" = H because H is real. The third 
line is a consequence of H being Hermitian. □ 

The quantum XYZ Hamiltonian has time reversal symmetry with respect to the spin-flip 15. We next demon- 
strate how to build up entanglement with such a system. Consider a collection of n qubits laid out in a cyclic array 
interacting under the Ising class of Hamiltonians given by Hxyz with J x = J y = 0: Hj s = L/=i ^;°~/' CJ }+i' wriere 
we identify o : n+l = a\ . 

The eigenvalues are given by {Xj} = {J z (n — 2Y,k jk © jk+\)\ j = hh ■ ■ ■ jn} (L61), where the addition is done 
modulo 2 over the components jt of the binary expansion of j. For n even, each eigenvalue Xj is paired with 
another of opposite sign and in particular, Xq = —Xn-i with \Xq\ = n\J z \ = X max . The concurrence spectrum of 
u = e~' H,s ' is composed of complex conjugate pairs and the concurrence capacity k„(m) may be computed explic- 
itly. Then k„(m) = max {iL^To 1 a 2 j e ~ 2 ' X ''\ > z'z — l,z T z — 0} where z — T!J=o a j\J)> P er Equation 14. Maximum 
capacity is obtained when the convex hull condition is satisfied which occurs precisely when the concurrence 
spectrum extends outside the right half of the complex plane. The minimum time at which this occurs is given by 

e -2in, in = ■ or tmin = ^4^1 = n/4n \ Jz i 

The existence of a time reversal symmetry in the interaction between qubits gives us important information 
about the nature of quantum correlations in the energy eigenstates. Applying Theorem V.7, we immediately find 
that the ground state of a Hamiltonian H with time reversal symmetry has maximum n-concurrence if it is unique. 
Examples of interactions satisfying these conditions are the XYZ Hamiltonian with (J x = /,. = /-=/> 0), denoted 
the XXX Hamiltonian, and the XY Hamiltonian (J x = J y ,J z = 0) (L61). In particular, the XXX Hamiltonian with 
/ > has been shown to have non-degenerate ground states in any number of dimensions, with or without periodic 
boundary conditions, provided the underlying lattice has a reflection symmetry about some plane (ibid.) 
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To illustrate this phenomenon we consider what happens when the time reversal symmetry is broken by adding 
a time-antisymmetric term to the XY Hamiltonian: 

h = ptpw+i + V o ^+' } + H (4o) 

where o" +1 = of. The presence of the linear term proportional to the total spin projection operator S z = YJj=\ 
breaks the time reversal symmetry so that iH ^ p when h z ^ 0. For zero magnetic field and < g < 1, the Hamil- 
tonian is time reversal symmetric and the ground state is non-degenerate meaning the concurrence is maximal. In 
the isotropic case (g = 0), the Hamiltonian commutes with S z and eigenstates are independent of h z . For magnetic 
field strengths below some critical value, \h\ < h crlt the ground state corresponds to an eigenstate with eigenvalue 
s z = of the operator S z . This ground state has maximal concurrence. For \h z \ > h cr u, the ground state corresponds 
to an eigenvalue s z ^ and the concurrence is zero (BB04a). 

VI. CONCLUSIONS 

We show that the odd-qubit concurrence canonical decomposition admits generalizations of all constructions 
studied on the even qubit CCD. In particular, a generalized pairwise concurrence capacity may be defined, and 
the operators for which this is maximal are characterized by a convex hull condition on the concurrence spec- 
trum. Again for an odd number of qubits, we find that for large odd n most unitaries have maximal concurrence 
capacities. Moreover, we provide an explicit algorithm for computing the odd-qubit CCD. 

These advances are complemented by new interpretation of the original inputs to the G = KAK theorem which 
define the CCD. Specifically, they may be rewritten in terms of time reversal symmetry 13 which is the spin- 
flip in « quantum bits, and the CCD is best understood in terms of such symmetries. For example, the odd- 
qubit CCD is a type AH KAK decomposition, and as such much have degenerate eigenvalues. In fact, this 
recaptures Kramers' degeneracy for the odd-qubit spin-flip, and a more careful study of the arguments reveals a 
Kramers' nondegeneracy: Nondegenerate eigenstates of 13 time reversal symmetric Hamiltonians only exist when 
the number of quantum bits is even and moreover must be highly entangled. Specifically, such \X) are highly 
entangled in the sense that the concurrence C„(|A,)) = | (A,|LS|A,) | = 1. Finally, the polar decomposition extracted 
from the CCD in the usual way accomplishes the following: any unitary n-qubit evolution is a product of precisely 
one time reversal symmetric and one time reversal antisymmetric evolution. 
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